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Abstract 

We numerically investigate the influence of self-attraction on the critical behaviour 
of a polymer in two dimensions, by means of an analysis of finite-size results of 
transfer-matrix calculations. The transfer matrix is constructed on the basis of the 
0(n) loop model in the limit n — > 0. It yields finite-size results for the magnetic 
correlation length of systems with a cylindrical geometry. A comparison with the 
predictions of finite-size scaling enables us to obtain information about the phase di- 
agram as a function of the chemical potential of the loop segments and the strength 
of the attractive potential. Results for the magnetic scaling dimension can be in- 
terpreted in terms of known universality classes. In particular, when the attractive 
potential is increased, we observe the crossover between polymer critical behaviour 
of the self-avoiding walk type to behaviour described earlier for the theta point. 



1 Introduction 



A useful formulation of the polymer problem can be given in terms of the 
0{n) loop model on a lattice, in the limit — * 0. In this formulation, loops 
have a weight n so that the number of loops is minimized as n — 0. In the 
simplest case, only one other parameter plays a role: a weight factor a may 
be assigned to each loop segment, in order to control the density of the loop 
configurations. A loop segment is the part of a loop that covers precisely one 
lattice edge. Loops do not intersect. For small a the vacuum state - without 
any loops - is stable; in contrast, for large a, most of the lattice edges are 
covered by a loop. The critical point separating the vacuum state and the 
dense phase has already been explored in detail; it is found to belong to the 
same universality class as the SAW (self-avoiding walk). 
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Exact results for the 0(n) loop model on the honeycomb lattice have been 
obtained [1-4] that reveal two branches of critical points (two points for each n 
in the range between and 2). One of these branches is interpreted as a critical 
point separating the disordered and long-range phases; the other branch is 
supposed to describe the ordered phase. This branch has algebraically decaying 
correlations for general n; in this sense it still qualifies as a critical branch. 

This picture applies to the case where the loop segments do not interact. A 
natural extension of the noninteracting loop model is to adjust the vertex 
weights in a way representing attractive forces between loop segments. Such 
interactions, if they arc of a sufficient strength, will influence the character of 
the phase transition between the vacuum and the dense phase. In the case of 
strong attractive forces, the critical state, which is relatively dilute, becomes 
unstable, so that the phase transition becomes first-order. 

The higher order critical point, separating the first-order transition from the 
continuous one, describes a polymer on the verge of collapse, and is called the 
theta point. An exactly known critical point of an n = loop model with 
vacancies on the honeycomb lattice, described by Duplantier and Saleur [5], 
belongs to the universality class of this theta transition. 

The formulation of the square 0{n) loop model given in Ref. [6] contains three 
adjustable parameters. One of these applies to the weight of vertices visited 
twice by a loop. There, both loop parts make 90° bends such that they collide 
but do not intersect. Thus, by varying this vertex weight one can tune the 
attractive forces. 

For this three-parameter model, there arc several branches of exactly known 
critical points [6-9]. Two of these have the same universal properties as those 
mentioned above for the honeycomb lattice. Two other branches are associated 
with a combination of Ising-like and 0(n) critical behaviour. The fifth branch, 
which was dubbed 'branch 0', is related with the dense phase of the 0(n + 1) 
model [6] and contains an intcgrable critical point [6,10] similar to the above- 
mentioned theta point [5]. In particular, the 0(n = 0) point has been shown 
to be equivalent to interacting walks on the Manhattan lattice at the theta 
point [11], with a set of exponents in agreement with those proposed for the 
theta transition [5], if one associates the magnetic dimension Xm = rj/2 — 
[5] with the critical dimension x^t,! which was discussed in Ref. [6]. 

Thus, the subset of the parameter space contains a considerable variety of 
critical points for which exact information is available [6-12]. Nevertheless, it 
covers only a small part of the parameter space. Therefore, it is of interest 
to explore the critical surface by numerical means. In particular, here we 
investigate the influence of attractive forces between the loop segments in the 
square 0(n) loop model, in the hope to observe the crossover between the 
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normal 0{n — 0) or SAW-like critical point and the proposed theta point. 

In Section 2 we briefly introduce the 0(n) loop model, and the transfer ma- 
trix used for the numerical calculations. Section 3 contains an analysis and 
a discussion of the numerical work. For this purpose, it also summarizes the 
flnite-size scaling formulas and an element of the theory of conformal invari- 
ance used in the analysis. 



2 The model and its transfer matrix 

The loop representation of the 0(n) model yields an expression for the par- 
tition function in terms of a sum over all graphs consisting of closed loops, 
covering a subset of the lattice edges [13,1]. The simplest case, with only two 
parameters, the loop weight n and the weight factor a per bond, has been 
investigated for the 0(n) loop model on the honeycomb lattice. In addition 
to exact results, conclusions based on a numerical investigation are available. 
Finite-size scaling of transfer-matrix results can be applied in parts of the pa- 
rameter space where no exact information is known. These results [14] have 
confirmed the interpretation of the two critical branches as mentioned above. 

In this work, we apply a similar method to the square lattice. The square 
0(n) loop model is now formulated in terms of loop weights n, bond weights 
a, and an additional weight factor p associated with each straight vertex, and 
an additional weight factor q associated with each 'collision'. Thus, the vertex 

weights are wq = 1 for an empty vertex (not visited by a loop), Wi = a for a 
vertex visited once by a loop making a 90° bend, W2 = ap for a vertex visited 
once by a straight loop part, and W3 = a^q for a vertex visited twice (see Fig. 
1). Using these weights, the partition function becomes 

Z = Y,<'w^'w^'n^\ (1) 
e 

where the sum is over all loop configurations The graph Q consists of Ni 
nonintersecting loops, with Ni {i — 1,2, 3) the total number of vertices of the 
type indicated by index i. 

A transfer matrix can be constructed for this loop model; a detailed descrip- 
tion is given in Ref . [6] . The transfer matrix can be seen as an operation that 
builds up a cylinder on which the lattice 0(n) model is wrapped, by adding 
one circular row, and thus increasing the length of the cylinder by one unit. 
The presence of the open end of the cylinder allows non-empty bond configu- 
rations even for n = 0, where closed loops are actually excluded. The transfer 
matrix indices are a numerical coding of 'connectivities': the way in which the 
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1 a a ap a^q 

Fig. 1. Vertex weights for the 0(n) loop model on the square lattice. In the present 
work we choose p = 0, so that 180° vertices are excluded. 

dangling bonds, or loop segments, are connected at the end of the cylinder. 
The allowed connectivities are 'well nested': no four occupied edges can be 
crosswise connected. This property is a consequence of the absence of inter- 
sections, and greatly restricts the number of connectivities. A sparse-matrix 
decomposition [6] allowed us to obtain transfer-matrix eigenvalues for systems 
up to linear size L — 12 using only modest computational resources. 

These transfer matrix eigenvalues are meaningful because they determine the 
free energy and the length scales determining the decay of the correlation 
functions along the cylinder. In the general case, the free energy per unit of 
area / of a model on an infinitely long cylinder with finite size L is determined 

by 

/(L) = L-MnAi°\ (2) 

where A^'* is the largest eigenvalue of the transfer matrix. In the paramagnetic 
phase of the 0{n — 0) loop model, as well as on the critical line separating 
the ordered phase, A^°^ = 1. The corresponding eigenvector is dominated by 
the vacuum. This eigenvalue 1 persists in the ordered phase, but there it is 
no longer the largest eigenvalue. For reasons of continuity, we maintain the 
notation A^"* = 1. 

For practical reasons, the set of connectivities is split into two disjoint subsets: 
the even subset, where all the dangling bonds are connected pairwise; and the 
odd subset where, in addition, a single, unpaired dangling bond occurs. The 
vacuum-dominated leading eigenvector naturally occurs in the even sector. 
The odd sector contains a line of covered bonds, running along the length 
direction of the cylinder. Such a hne is not a part of a closed loop, and does 
not occur in the graphs Q of Eq. (1). Therefore, the odd sector does not 
contribute to the partition sum. 

However, the odd connectivities are important for the calculation of the mag- 
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netic correlation length of the 0(n) model. The graphs containing, in addition 
with closed loops, a single line of covered bonds connecting two points are pre- 
cisely those describing the magnetic correlation function between 0(n) spins 
located on those points. It is expressed as Z' jZ^ where Z is given by Eq. (1), 
and Z' by the same equation but with Q replaced by Q' representing all graphs 
consisting of closed loops plus a single line connecting the two spins. 

Thus, the inverse magnetic correlation length ^~^(L) can be expressed in terms 
of A^^-* and A^^-* which denotes the largest eigenvalue in the odd sector: 

r^(i>)=ln(Af/AW). (3) 

A convenient quantity in the finite-size analysis is the 'scaled gap' x(it, V) — 
L/\2n^{u, L)] whose arguments include, besides the finite size L, also a temper 
ature-like parameter u. For instance we may choose u = a—a^ where denotes 
the critical value of the bond weight a. But u may also be chosen to parametrize 
the distance to fixed points located on the critical surface. In the vicinity of a 
renormalization fixed point, finite-size scaling leads to the equation 

x{u, L)^Xh + ^Ly-u[dC\u, l)/du]u=o + ..., (4) 

ZTT 

where Xh is the magnetic scaling dimension [15] of the corresponding fixed 
point, and i/u the renormalization exponent associated with u, governing the 
flow to or away from the fixed point. Thus, if i/^ < (> 0), x{u, L) converges to 
(diverges from) xy^ with increasing L. This allows us to analyse the behaviour 
of the finite-size results x{u, L) in the light of the phase diagram. 



3 Results 

Before starting the actual transfer- matrix calculations at n = 0, we summarize 
the role of the three parameters. First, the bond weight a (or wi) is adjusted 
in order to find the critical point. Expressing the scaled gap in a, ignoring the 
irrelevant scaling fields, and using finite-size scaling of the correlation length, 
the critical point Oc is determined by numerical data for two subsequent even 
finite sizes L and L -\-2: 

x{ac, L) = x{ac, L + 2). (5) 

Only even sizes are used because, in general, the odd systems display differ- 
ent finite-size amplitudes [6] . Corrections to scaling introduce deviations from 
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Eq. (5) so that extrapolation of the finite-size estimates of Oc was performed. 
For details, see e.g. Ref. [16]. 

The parameter p is here important with regard to the Ising-like degrees of 
freedom that play a role in the square loop model [6]. In order to clarify 
this point, we introduce an Ising spin on each elementary face, such that two 
neighboring spins are equal only when a loop passes between them. Clearly 
when p = a loop is adjacent to spins of one sign only, signalling a broken 
Ising symmetry. This means that an Ising-like ordering transition occurs when 
p is lowered. 

In the present work, we wish to focus on the collapsing polymer problem, 
while avoiding the interfering effects associated with Ising-hke ordering. A 
possible way to circumvent these effects is to exclude type 2 vertices by putting 
p = W2 = 0, so that the Ising degrees of freedom are already frozen even in the 
relatively dilute SAW-like critical state. An increasing polymer density does 
not further affect this Ising ordering. 

Thus, we scanned the parameter space using wi = a, W2 = 0, and = qa^ 
where the parameter q governs the attraction between the loop segments. This 
choice of parameters includes the integrable theta point mentioned above, for 
which a = 0.5, q = 2. 

Using finite size parameters L = 2, 4, 6, 8 and 10, Eq. (5) was solved for q = 0, 
1.0, 1.5, 1.8, 1.9, 2.0, 2.1, 2.2, 2.5, 3.0 and 4.0. The extrapolated critical points 
are shown in Table 1. 

The value ac = 0.63860(5) at g = gives the connective constant 1/a,. ~ 
1.565... for SAWs with a 90° turn at every step, in good agreement with 
the estimate 1.5657(19) [17]. The scaled gaps at the intersections are plotted 
versus the system size L in Fig. 2. These data are to be compared to Eq. (4) 
where we may interpret i as a field parametrizing the critical line; the leading 
temperature field vanishes in effect for the solution of Eq. (5). The field along 
the critical line is expected to be irrelevant in the vicinity of the SAW-like 
fixed point; at a higher critical point, it is expected to be relevant. Indeed, 
ior q < 2 we observe, for increasing L, a converging trend towards the exactly 
known value Xh — 5/48 = 0.104166 ... for the SAW model. This convergence 
reflects the stabihty of the SAW-hke flxed point. For q = 2 the finite-size 
data converge well to Xh — 1/4, in agreement with the known value at the 
integrable point. 

For q > 2 the phase transition becomes discontinuous as a function of the 

temperature-like parameter a, as is revealed by an intersection of the two 
largest eigenvalues Aq = 1 and A2 of the transfer matrix in the even sector. 
For the interpretation of the results shown in Fig. 2 it is important to note that 
we used Ao = 1 in Eq. (3), even where other eigenvalues exceed 1; this occurs 
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Table 1 

Numerical results for the critical value of the bond weight a for 0(n = 0) loop 
models with different values of the loop-loop interaction parameter q (see column 
1). For q < 1.5 we assumed a correction term proportional to L~^/^, for q > 1.5 
one proportional to L~^. The rightmost column shows the estimated value for the 
magnetic dimension Xh where extrapolation was possible. The extrapolations for 
q = 2 are in a good agreement with the exact values Oc = 1/2 and Xh = 1/4. 



q 






0.0 


0.63860 


(5) 


0.1045 


(5) 


1.0 


0.5769 


(1) 


0.104 


(1) 


1.5 


0.5399 


(1) 


0.11 


(1) 


1.8 


0.5165 


(1) 






1.9 


0.5084 


(1) 






2.0 


0.5001 


(1) 


0.2500 


(2) 


2.1 


0.4917 


(1) 






2.2 


0.4833 


(2) 






2.5 


0.4575 


(5) 






3.0 


0.421 


(1) 






4.0 


0.368 


(1) 







for q > 2 even at the line of phase transitions. Using this analytic continuation 
we avoid irregularities caused by the intersections mentioned. The increasing 
trend of the data in Fig. 2 for g > 2 reflects the instability of the theta-hke 
fixed point. 

However, the lines connecting the finite-size results in Fig. 2 are running almost 
horizontally for g f« 2; the finite-size dependence of the scaled gap is rather 
weak. This corresponds with a small value of the exponent in Eq. (4) when 
applied to the theta point. Indeed, from Coulomb gas arguments [5,19] one 
expects Uu — 3/4. 

Also in the vicinity of the stable SAW-like fixed point, the finite-size depen- 
dence of the data shown in Fig. 2 is rather weak. This stands in a strong 
contrast with the rapid convergence observed [18] for 0(n = 0) loop models 
with p — 1, which can be interpreted [16] in terms of an irrelevant exponent 
Ui = —2. The present slow trends near the SAW-like fixed points may be ex- 
plained by an irrelevant exponent —11/12, in analogy with the case of 'trails' 
where loops are allowed to intersect [16]. In the Coulomb gas representation, 
loop intersections and collisions correspond with the same diagrams. 
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Fig. 2. Finite-size dependence of the scaled gap of 0(n) models. These data are 
taken at values of the bond weight a solving the scaling equation for the correlation 
length for two subsequent system sizes. The lines connecting the data points are for 
visual aid only. Starting from below, the lines apply to g = 0, 1.0, 1.5, 1.8, 1.9, 2.0 
(dashed line), 2.1, 2.2, 2.5, 3.0 and 4.0 respectively. 

We have performed extrapolations on the basis of the numerical solutions of 
Eq. (5). According to finite-size scaling, the sohitions a{q,L) obtained from 
finite sizes L and L ~ 2 behave as a{q, L) = a{q) + bL^^'^^ + • • ■ where yt is 
the leading temperature exponent; yt = 4/3 at the SAW-like fixed point, and 
yt = 7/4 at the theta point. For the second temperature-like exponent we may 
expect yu — —11/12 and y„ = 3/4 respectively. On the basis of our limited 
range of finite sizes, it was not possible to well determine y„ independently 
for all q. Only in the case g = we could confirm that yu — Ut ~ —2.2 as 
expected. For q = 2.1 and 2.2 we find results for y^ — yt consistent with — 1 
but the data show that further corrections are important. Thus we assumed 
yu~yt — —9/4 for g < 1.5 and — = — 1 for g > 1.5, and extrapolated the 
finite-size estimates of the critical points accordingly. The results are shown 
in Table 1. More details of the fitting procedure can be found in Ref. [16]. 

Since the solutions of Eq. (5) converge to the critical point, the scaled gaps 
converge to the corresponding magnetic scaling dimension. For those cases 
where the finite-size convergence was sufficiently fast, extrapolated results are 
included in Table 1. The results for g = 0, 1 and 1.5 are in a good agreement 
with the expected 0(n = 0) magnetic dimension Xh = 5/48 = 0.104166... 
Thus, the fact that the Ising degree of freedom is frozen docs not seem to 
alter the universal properties. For q = 2 our result for the magnetic dimension 



8 



agrees well with xu = 1/4 as found earlier for this integrable point of the square 
0(n = 0) model [6,10]. For larger values of q the finite-size data assume a 
monotonically increasing trend with L, in accordance with earlier observations 
at first-order transitions, for instance that of the Potts model with a large 
number of states [20] . 

The phase diagram of the model, summarizing our findings, is shown in Fig. 3. 
The vacuum state occurs near the origin, and a dense phase at large values of 
Wi and W3. These phases are separated by a line of phase transitions of which 
the upper part is first order, and the right hand part is critical, with the theta 
point separating both parts. 




Fig. 3. Phase diagram of the 0(n = 0) model for W2 = 0. The Une of phase transi- 
tions has a first-order part (soUd line) and a continuous part (broken line) . This line 
is an interpolation between the numerical results (open circles), the exactly known 
theta point at = 1/2 (black circle) and the exactly known point at wi = (black 
square) . 



In conclusion, our data confirm that the introduction of sufficiently strong 
attractive forces in the critical 0(n = 0) model leads to the collapse of a 
polymer in two dimensions. The critical exponents at the point of collapse 
are in agreement with the predictions of Duplantier and Saleur [5] provided 
we avoid the onset of Ising-like ordering which is known to lead to a different 
universality class [6]. 
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